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Introduction 

The insolation parameters of the Earth depend on its orbital parameters and of the pre- 
cession and obliquity. Until 1988, the usually adopted solution for paleoclimate computation 
consisted in (Bretagnon, 1974) for the orbital elements of the Earth, which was completed by 
(Berger, 1976) for the computation of the precession and obliquity of the Earth. In 1988, I 
issued a solution for the orbital elements of the Earth, w^hich was obtained in a new manner, 
gathering huge analytical computations and numerical integration (Laskar, 198S). In this 
solution, w T hich will be denoted La88, the precession and obliquity quantities necessary for 
paleoclimate computations were integrated at the same time, w r hich insure good consistancy 
of the solutions. Unfortunately, due to various factors, this latter solution for the precession 
and obliquity w T as not wndely distributed (Berger, Loutre, Laskar, 19S8). On the other side, 
the orbital part of the solution La88 for the Earth, w^as used in (Berger and Loutre, 1991) 
to derive an other solution for precession and obliquity, aimed to climate computations. I 
also issued a new solution (La90) which presents some slight improvements with respect to 
the previous one (Laskar, 1990). As previously, this solution contains orbital, precessional 
and obliquity variables. In order to make it wddely available, it was distributed during this 
meeting on magnetic support and can also be obtain directly by request to the author at 
laskar®friap51 .bitnet. In the present talk, I will discuss the main features of this new 
solution. 




The orbital solution La90 


The orbital solution is obtained by the numerical integration of an extended averaged system, 
which represents the mean evolution of the orbits of the planets. All the S main planets of the 
solar system are taken into account, as well as lunar and relativistic main perturbations. The 
use of numerical integration for the computation of the solution of the secular system is one of 
the reasons for the good quality of this solution, which can be checked by comparison with the 
available ephemeris over a short time scale (Laskar, 1985, 1986, 1988). In (Laskar, 1988), the 
solution LaSS was represented in quasi-periodic form over 10 Myr, but these representations 
are slowly convergent, which prevents a good accuracy for the solution. Later on, I understood 
fully the reason of this slow convergence, which is due to the presence of multiple resonances 
in the secular system of the inner solar system. Due to these resonances, the motion of the 
solar system is chaotic, and not quasi-periodic, as was demonstrated by the computation of 
its Lyapunov exponents which reachs 1/(5 Myr) (Laskar, 19S9). This implies that it is not 
possible to give any precise solution for the motion of the Earth over more than about 100 
Myr, and most probably, ephemeris can only be given for about 10 Myr with good precision. 
Several integrations of the secular system of the solar system were made over 200 Myr and 
400 Myr. The origin of the chaotic behaviour was identified, and is due to the presence of 
secular resonances in the inner solar system (Laskar, 1990). With a new numerical method, 
it was possible to show' that the chaotic zones w-ere the solar system belongs is large in the 
directions of the proper modes related to the inner planets. This is an indication that these 
results are stable against small changes of initial conditions or model. 


One of the main consequence of interest for paleoclimate studies, is the fact that the main 
frequencies of the orbital motion of the planets can no longer be considered as constants, but 
are slowly evolving with time. The measured shift in frequency amount about 0.2 arcsec/vear 
over 200 Myr for g 3 and g 4 w'hile for g 5 it is of only 0.00002 arcsec/year. It should be stressed 
that, as the motion is chaotic, the computed solution cannot be considered as close to the 
real solution over more than a few 10 Myr, but it is reasonable to assume that the drift 
m frequency is of the same order of magnitude over the 200 Myr. One should nevertheless 


mention that over this time span, the change of the precession main frequency, due to tidal 
effects in the Earth-Moon system, are probably more important (Berger et ai, 1989). 

Since, direct numerical integrations over 3Myr backward and forw^ard w r ere issued by (Quinn 
et ai, 1991), including also solutions for precession and obliquity (QTD6). The orbital 
solutions have been compared with La90. The tw'o solutions present very small discrepancies 
over 3 Myr, and the existence of the secular resonances in the inner solar system is confirmed 
(Laskar, Quinn, Tremaine, 1991). The very close agreement of the two orbital solutions over 
3Myr gives the confirmation that the Earth parameters are now very well known over this 
time span, and insure that the orbital solution La90 can be used w r ith confidence over 10 Myr 
for paleoclimate use. The precession and obliquity solutions present some small discrepancies 

which are probably due to the presence of the tidal effect in the Earth-Moon system in the 
QTD6 solution. 
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Precession and Obliquity in the La88 and La90 solutions 


The precession quantities are completely determined by the two motions of the equatorial 
and ecliptic pole. The motion of the ecliptic is given by the secular theory La90 (Laskar 
1990); the precession quantities are integrated at the same time, using the equations of the 
theory of the rigid Earth of Kinoshita (Kinoshita 1975, 1977, Laskar, 1986). The equations 
for the general precession in longitude pa , and for the obliquity of the date e are then 


with 
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where p = sin(z/2) sin(fl), q = sin(t/2) cos(fl), (i is the inclination of the Earth with respect 
to a fixed ecliptic, and f l the longitude of the node). J?(e) is the secular term due to the 
direct lunisolar perturbations. The quantities M 0 , Mi, M 2 , M 3 , So, and S 2 depend only on the 
orbital elements of the Moon and the Sun. The pincipal moments of inertia of the Earth are 
denoted by A, B , and C, and the angular velocity of the Earth is u . The masses of the Sun, 
the Earth, and the Moon are denoted by m©, mg, and rriM ; the sideral motion of the Sun and 
of the Moon by n© and um ; and the mean motion of the node of the Moon by uq. The other 
terms present in equation (23) represent the effects of the secular variation of the ecliptic, 
caused by the secular planetary perturbations. The numerical values of Mo, Mi, M 2 , M3 are 
given in (Kinoshita, 1977): 

M 0 = 496303.3 x 10" 6 
Mi = -20.7 x 10" 6 
M 2 = -0.1 x 10~ 6 
M 3 = 3020.2 x 10~ 6 
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and from (Laskar, 1986), 


So = ^(1 - e 2 )" 3/2 - 0.422 x 10 -6 

The following numerical values, were also used (see Laskar, 1986 for complete references) : 

474 659981.59757 "/yr 
17325593 .4318 "/yr 
-69 679.193 6222 "/yr 
384 747980.645 m 
0.017202098 95 
328 900.5 
332 946.0 

The quantity p g is the geodesic precession due to the general relativity, 

p g = 0.019188 ; 7yr 

The value of the dynamical ellipticity Ed = (2 C — A — B)/2C = 0.0032S005 is obtained by 
adjustment at the origine J2000 to the values of the speed of precession and obliquity give 
by the IAU (Grenoble, 1976) : 

p = 50.290966'Vyr 
e 0 = 23°26'21"448 

For t = 0 , we have Pa = 0 ,£ = eo , i = f2 = 0 ,and thus : 


u = 

tim = 

rcn = 

dM = 
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mQ/(m E + rn M ) = 
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= R(e o ) — 2 p t=0 - p g cot e 0 

(=0 

These formula for precession gives a solution for precession and obliquity in agreement with 
the requirements of high accurate ephemeris (Laskar, 1986) and are thus best suited for 
paleoclimate computations. 


dpA 

dt 


Descripion of the files of the solution La90 


The different files which are distributed are 
ear0m5.dat orbital elements 0 to -5 Myr 

pre0m5.dat obliquity and precession 0 to -5 Myr 

ear5ml0.dat orbital elements -5 to -10 Myr 

pre5ml0.dat obliquity and precession -5 to -10 Myr 

All the elements are referred to the ecliptic and equinox J2000. The starting date is J2000. 
All angles are in radians 

The ear0m5.dat and ear5ml0.dat files contains T, k,h,q, p 
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T = time from J2000 in lOOOyr 
k = e cos(ro) 
h = e sin(ro) 
q = sin(i/2) cos(fl) 
p = sin(i/2) sin(fl) 

where e is the eccentricity of the Earth, i denotes the inclination of the Earth, zzj the longitude 
of perihelion, and S7 the longitude of the node of the Earth with respect to the fixed ecliptic 
and equinox J2000. 

The pre0m5.dat and pre5ml0.dat files contains T,e,p A where 

T is the time from J2000 in lOOOyr, £ is the obliquity of the date (mean equator of date with 
respect to mean ecliptic of date), and p A is the general precession in longitude. 

All quantities generally used for climate studies are derived easily from these fundamental 
quatities. The eccentricity of the Earth is obtained by 


e = y/k^ + h 2 

And the longitude of perihelion from moving equinox is 


UJ* = G7 + PA 

The climatic precession e cos(ui*) is thus equal to the real part of z*exp(ip A ), where z = k+?h. 
Conclusions 

The present solution La90 for orbital and precession elements for the Earth over can be used 
for as an ephemeris paleclimate computations over 10 Myr. On this time span, the fact that 
the motion of the solar system is chaotic is not perceptible on this level of precision. This 
new solution present some improvements from my previous solution LaS8 (Laskar 19S8). The 
orbital solution is in very good agreement with the recent numerical integrations of (Quinn et 
al . , 1991) . Over the geological time scales exceeding 100 Myr, there is no hope to obtain such 
an ephemeris, due to the chaotic behaviour of the solar system, but one can assume that for a 
few 100 Myr, the slow diffusion of the astronomical frequencies observed during the 200 Myr 
integrations remains of about the same order. This should not be granted for billion years 
computations and more computations on the diffusion process in the solar system are clearly 
needed. More, due to the changes in frequencies, and presence of secular resonances in the 
orbital forcing, more extended studies of the resulting influence on the rotational evolution of 
the Earth should be done. Studies on the long term evolution of the Earth-Moon system are 
also needed in order to improve our knowledge of the rotational evolution of the Earth. The 
geological records, assuming that a good accuracy in the determination of the fundamental 
frequencies could be achieved, would then be the only possible observations for tracking the 
long term past evolution of the solar system for time span longer than 100 Myr. 
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